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The brane world model proposed by Dvali-Gabadadze-Porrati realizes self-acceleration of the 
universe. However, it is known that this cosmological solution contains a spin-2 ghost. We study 
the possibility of avoiding the appearance of the ghost by slightly modifying the model via the 
introduction of a second brane. First, we consider a simple model without stabilization of the brane 
separation. By changing the separation between the branes, we find that we can erase the spin-2 
ghost. However, this can be done only at the expense of the appearance of a spin-0 ghost instead. 
We discuss why these two different types of ghosts are correlated. Then, we examine a model with 
stabilization of the brane separation. Even in this case, we find that the correlation between spin-0 
and spin-2 ghosts remains. As a result we find that we cannot avoid the appearance of a ghost by 
introducing a second brane in the model. 

o : 

I. INTRODUCTION 

x> : 

D . n, 

The observed present-day accelerated expansion of the universe [1| is one of the hottest topics in cosmology. There 

are two aspects in this issue: The first is to explain the late-time accelerated expansion of the universe, and the second 

is to explain why the accelerated expansion began only recently. There have been many attempts to modify standard 

cosmological models in order to address these problems. As far as we know, none of these give a natural solution to 

(f-) \ the second problem. To give a solution to the second problem, one may need to combine ideas which solve the first 
■ problem with the anthropic principle [3| ■ 

C\| ' In this paper, our interest is in the first problem. Most of ideas that solve the first problem are modifications of the 
scalar (spin-0) sector of the cosmological model, such as introduction of a cosmological constant Q or quintessence [|J . 
But there is another direction which has not been explored much so far. That is modifying the gravity theory in the 

~j , spin-2 sector. 

The simplest model with a modified spin-2 sector involves massive gravity [5|. Since models which have terms 
quadratic in metric perturbations can be regarded as a massive gravity theories, most modified gravity models fall 
""^ • into this category. If we introduce the mass of the graviton to explain the accelerated expansion of the universe, 
its value would be the same order as the present value of the Hubble parameter, m g ~ H . However, it is known 
JL , that a spin-2 graviton with mass in the range < m 2 < 2H 2 in de Sitter background has a ghost excitation in 
its helicity-0 component Q. The Hubble parameter is larger in an earlier epoch of the universe. Hence it seems 
difficult to avoid the appearance of a ghost throughout the evolution of the universe in simple models that attempt to 
L* ' explain the present-day accelerated expansion. Here we would like to raise a simple question: Can we build a model, 
. £h ! incorporating modification of the spin-2 secton that explains the accelerated cosmic acceleration and is ghost- free? 
^ . The DGP brane world model (See Ref. [| for a recent review paper), in which the 4D Einstein-Hilbert action 
is assumed to be induced on the brane, is a mechanism to realize the late-time accelerated expansion of the Universe 
without introducing additional matter Q. Such a solution is called the self-accelerating branch. This is in contrast 
to the normal branch, where there is no cosmic expansion in the limit of vanishing matter energy density. There 
are many investigations of the DGP model, such as spherically symmetric solutions [l(| QjJ, analyses of structure 
formation [12] and shock wave limits [l4| . In quantizing the theory, we face the issue of strong coupling [IH, [h| , which 
is the problem that nonlinear effects of the scalar gravitational interaction become important below some unexpectedly 
large length scale. Moreover, the existence of a ghost excitation has been pointed out in the self-accelerating branch 
of the DGP brane world model Oil 03 ■ We here discuss this ghost problem. By considering perturbations around 
this solution, a Kaluza-Klein tower of massive gravitons in the four dimensional effective theory is obtained. The 
lowest mass m satisfies m 2 = 2H 2 [H, If the analogy to the massive gravity theory holds, this would mean that 
the self-accelerating branch of the DGP model does not have a ghost excitation. However, it a detailed analysis has 
shown that there is still a ghost in this model [ItJ • 

In this paper, based on the DGP brane world model, we attempt to build a ghost-free model which simultaneously 
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explains the accelerated expansion of the universe due to effects in the spin-2 sector. The first idea is to increase 
the lowest graviton mass by introducing a boundary brane. We call this the two-brane model, which is successful in 
making all graviton masses satisfy 2H 2 < m 2 . However, a spin-0 excitation, which originates from the brane bending 
degrees of freedom, is transmuted into a ghost. We will find that a ghost excitation cannot be eliminated in this 
two-branes model. 

Usually spin-2 modes are decoupled from spin-0 modes at the level of linear perturbations. Therefore it seems 
mysterious that the disappearance of the spin-2 ghost coincides with the appearance of a spin-0 ghost. The magic 
is in that a spin-2 mode can be obtained from a spin-0 mode by applying a differential operator if and only if the 
mass squared takes a specific critical value. At the critical mass, a spin-2 mode can couple with a spin-0 mode. In 
addition, the critical mass corresponds to a threshold for the appearance of a ghost in the spin-2 sector. If the mass 
of a spin-2 mode is smaller than the critical mass, the spin-2 sector has a ghost excitation as previously mentioned. 
Moreover, a unique spin-0 mode that appears in the two-branes model has also the critical mass. Hence the spin-2 
mode can couple with the spin-0 mode when it crosses the critical mass. This fact partly explains why the ghost can 
be transferred between the spin-2 and spin-0 sectors. 

Then, we consider a model with stabilized brane separation by introducing a bulk scalar field into the two-branes 
model [H, [l9{ . Since the mass spectrum of spin-0 modes changes, the spin-0 mode at the critical mass does not 
survive in this stabilized two-branes model. Then, it is expected that the ghost will not be transferred between the 
spin-2 and spin-0 sectors. However, by studying this stabilized two-branes model in detail, we will find that this naive 
expectation is wrong. 

This paper is organized as follows. In £|TT]we discuss the two-branes model. In i jllll extension to the stabilized 
two-branes model is performed. We will summarize our results in £ )IVI 

II. TWO-BRANES MODEL 

As mentioned in Introduction, it is known that in de Sitter spacctime a massive spin-2 field whose mass squared 
is less than 2H 2 has a ghost [6(, where H is the Hubble parameter. The self-accelerating branch of the DGP brane 
world model @, Q has a spin-2 excitation whose mass squared is 2H 2 . In the usual massive graviton model with 
Fierz-Pauli mass term, this is a special case where a helicity-0 excitation disappears [f| |2(j. However, a detailed 
analysis reveals that there is a ghost mode in the self-accelerating branch of DGP model due to the existence of a 
spin-0 excitation [13] ■ Since the ghost arises due to a subtle mechanism in this model, it would be a natural question 
whether one can eliminate the ghost by a slight modification of this model or not. 

In this section we examine the two-branes model with a boundary brane in the bulk as one of the simplest extensions 
of original DGP scenario. Our initial expectation is as follows: The new boundary will shift the mass of spin-2 modes. 
Then we will be able to construct models with no spin-2 excitations with masses in the range that implies the existence 
of a ghost. 

However, we will find later that spin-0 excitations, which correspond to the brane bending degrees of freedom, carry 
a ghost mode instead. It will turn out to be impossible to eliminate the ghost by this simple two-branes extension. 
Below, we discuss how the ghost is transferred between the spin-2 and spin-0 sectors [22| . 

A. Background 

As mentioned above, in order to obtain a model that does not have spin-2 ghost modes, we introduce another brane 
into the DGP model. The action is 




where R, R± , g± , K±, t± and L m ± are the five dimensional Ricci scalar, four dimensional Ricci scalar, trace of 
the four dimensional induced metric g^ v , trace of the extrinsic curvature K llUl the brane-tension, and the matter 
Lagrangian on the (±)-branes, respectively. Assuming Zi symmetry across each brane, the junction conditions to be 
imposed at each brane are 



(2.2) 



with r c := k 2 /2n\. 
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We assume vacuum without any matter field, L ma = 0, as an unperturbed state. The five dimensional metric is 
given by 

ds 2 = dy 2 + a 2 (y)"f^dx^dx", (2.3) 

with two boundary four dimensional de Sitter (±)-branes, where 7^^ is the four dimensional de Sitter metric with 
unit curvature radius and the warp factor a(y) is given by 

a{y) = a+ + y. (2.4) 

The range of y-coordinate is from to (a_ — a + ). The positions of (+) and (— )-branes are at y — y + :— and 
at y = y_ := a_ — a+, respectively. In the present model, the value of the warp factor on the brane is related to 
the Hubble parameters evaluated on each brane H± as a± = 1/H±. Notice that the bulk geometry of the above 
solution is nothing but the five dimensional Minkowski spacetime written in the spherical Rindler coordinates. From 
the junction conditions (|2.2p . the Hubble parameters on the respective branes are related to the tension of the branes 
as 

±H± = r c Hl - —t ± , with H + > H_> 0. (2.5) 



B. Perturbation equations 



Now we study perturbations around the background mentioned above. Besides the transverse-traceless conditions, 
one can impose the conditions of vanishing {yy} and {?//i}-components of metric perturbations. Then the perturbed 
metric is given by 



ds 2 = dy 2 



+ (a 2 ltlv + h% T ^dx^dx u , 



with 



(2.6) 



(2.7) 



where V M is the covariant derivative operator associated with 7^ 
use ^ijiv The {/i^}-components of the Einstein equations in the bulk are 



a 2 l^u- 



xr ) -2e T) =-( a(4) -2)e/ ) , 



In raising or lowering Greek indices, we 

(2.8) 



where a prime " ' " denotes a partial differentiation with respect to y, and := 7 A " y V AI V„. 

In this gauge, branes do not in general reside at fixed values of y. We denote the locations of the branes by 
y = y± + £±\x ,i ). To study the boundary conditions, it is convenient to introduce two sets of new "Gaussian normal" 
coordinates associated with respective branes. Here "Gaussian normal" means coordinates such that {?/ /^-components 
vanish and the location of each brane is specified by a constant y-surface. Notice that the "Gaussian normal" 
coordinates associated with (+)-brane in general differ from those associated with (— )-brane. In these coordinates 
the junction conditions imposed on the branes is simple. We associate an over-bar "~" with the perturbation variables 
in these coordinates. 

The junction conditions for metric perturbations in the "Gaussian normal" coordinates are 



2rv 



{lis 



with 



X& :=-- 
^ 2 



□ (4) ^ - v^ch^i" - v v v a m a 



2%» (v a V^^j 



(2.9) 



(2.10) 



Here T( ± ) := j^T^ and X& is defined in the same manner. The generators of the gauge transformation from the 
"Gaussian normal" coordinates to the Newton gauge are given by 
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and 



tf ±) := x*-x*= -^ v d v e {±) {xt>)+e {±) {x»)- 
By this gauge transformation, the metric perturbations transform as 

h${y) = hZ T) (y) - 2a (1 - aH ± ) V„V^ ±) - 2a^£f ±) (x<>), (2.11) 

where we have chosen £^_js so that the second term in p. lip vanishes on the brane. Substituting (|2.11j) into the 
traceless part of Eq. 



the perturbed junction conditions for h^u imposed at y = y± becomes 
±hff* T 2H ± hZ T) + r c Hl (□«> - 2 ) = , 



with 



where 7+ := 7 All/ | y= j /+ . Combining the bulk equations and the junction conditions, we finally obtain 



(2.12) 



(2.13) 



l(TT) 



□ W- 2 



hZ T) iv) = E (-a^fc? - 2r c ^ (□«> - 2) <% - (2.14) 



<J=± 



where 



L (TT) ._ I* 4-1 



The trace part of Eq. (|2.9[l gives 



(□W + 4)^ : 



! ± 6iJ|(lT2r c ff±; 



T 



(±) 



(2.15) 



(2.16) 



The source Y,ffi should also satisfy the transverse conditions, which can be shown to be identical to the condition (|2.16|) 
by using the identity, 



which holds for an arbitrary scalar function Y. 



(2.17) 



C. Solution with source 



Equations (|2.14p are separable. We define a complete set of functions of y by (real-valued) eigen-functions of the 
eigenvalue equation associated with (|2.14p . 



^ l + 2r c £%-y CT ) 



<7=± 



i (y)=-L(^) Ui(y ) ) 



(2.18) 



where mf is the eigenvalue, which we find corresponds to cK 4 ) — 2 in comparison with Eq. (|2.14p . Then, thus defined 
eigen-functions are mutually orthogonal with respect to the inner product defined by 



{ui,Uj)^ := J> % [ l + 2r c £ 8{y - y a ) ) u t (y)u 3 (y), 



(2.19) 



where the integral is taken over a twofold covering coordinate y which ranges from — (y_ — y + ) to — y+) and 

is related to y by y = |y|. Both ends of the range of y are identified. By definition, the norm (iii, Ui)( TT ^ is positive 
definite. We normalize the eigen-functions so that they satisfy (ui,Uj)( TT ^ = 5{j. 
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l r T r F\ 

Expanding hpv in terms of these eigen-functions as 

^ T) -EM^K(y)' ( 2 - 2 °) 

3 

the equations of motion (12. 14[) reduce to 

i E ( d(4) ~ 2 ~ m f) + 2 ^ E *(» - ^)) *&V)«i(t/) = -^$5(y - y+), (2.21) 

where we have assumed that the source is only on the (+)-brane. Operating (^ dyui(y)x) on both sides of the above 
equations, we obtain 

(□W- 2 - m?)fc«(s") = -2 K 2 I]W^(y + )- (2-22) 
Hence, the solution for becomes 

i 1 

From (|2. 1 1|) with the aid of 12 . 16[) and (I2.23|) . the induced metric on the (+)-brane for the solution with source is 
given by 



U+) - 2/r 2 V tf(y+) T (+) I- T (+) 

i 1 \ 



3(^2)^-4^ J 



' 6ff 2 \i m ^ 2 / / D(4) 

where we have neglected terms which can be erased by a gauge transformation, and we have used the relation 

1 / 1 \ / 1 _f/n\ 1 



V M V„ - ~ 7/ ^ (4) ) Y = ( V„V„ - \ lllv U& ) -re ^ : - y, (2.25) 



□ ( 4 )-2-m 2 V 4 ,Ml/ y V 4 ,Ml/ yn(4)-m 2 + 6' 

which holds for an arbitrary scalar function Y. This relation is equivalent to the identity 



(v^V„nW - nWv.V,) y = -8 fv M V„ - ^7^ n(4) ) Y- (2-26) 



Note that we have also used the relation 

1 



(□(4) _2-m 2 )(n( 4 ) -4) m?-2 \aW -2-mf n( 4 ) 



(2.27) 



which cannot be applied when mf — 2. 

We stress that when mf = 2, there is no pathology present in the model, in contrast to the Fierz-Pauli case. In 
the Fierz-Pauli model, we cannot put the matter with T ^ if the squared mass is set to 2H 2 , because the equation 
of motion implies that T = (see Eq. (|A7|) in the Appendix). However, in the DGP brane world model there is no 
such a restriction. When the DGP brane world action is reduced to the form of Eq. (|2.24[) . there appears to be a 
similar pathology. However, this naive guess is wrong since Eq. (|2.24|) is derived by using Eq. (|2.27p . which is not 
correct when mf = 2. When we go back to more basic equations (|2.23[) and (|2.16p . it is easy to find that the matter 
with T ^ can be consistently introduced into the DGP brane world even if mf = 2. The seemingly pathological 
expression (|2.24p at mf — 2 is due to the decomposition to spin-2 and spin-0 perturbations (|2.27p . 



D. Effective action 

From the expression (12.24[) . we find that the physical gravitational degrees of freedom are composed of a Kaluza- 
Klein tower for the transverse-traceless part and a single mode corresponding to the trace part in K^} . We denote 

- (z) 

them by h^l and s, respectively. The induced metric perturbation is given by 

i 
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The kinetic term of the effective action is written as 



Ski n = [ d i Xy/^h^(n-2-rn 2 )h% 

+13 [ d A xy/^%s(n + 4) S: (2.28) 



where the indices are raised or lowered by using 7+ On the other hand, the action of the matter localized on the 
(+)-brane, expanded in powers of hpw, is written as 

Setter = (0" th OTder of h) + E \ J faV^Tffitfl 

i 

+ I /"d 4 xy=^T(+) S + 0(h 2 ). (2.29) 

Comparing the equations of motion derived from the variation of Skin + S ma tter with Eq. (|2.24[) . we can fix the 
coefficients, on and j3. Here a little subtlety arises when we take variation with respect to h^l, since this variable 
is constrained to be transverse-traceless. To avoid complication, here we simply consider traceless source; T^ + ' = 0, 
deferring detailed discussion to Appendix. Then, we find 



^ = Q 2 if V ( 2 - 30 ) 



As for the trace part, by comparing simply the poles at n 1 - 4 ) + 4 = 0, we obtain 

ml 



.M(j_ ( ^ + _ ir . +2 ( ? *i))-. 



The remaining tasks are to find the mass spectrum of the eigenvalue equation (|2.18[) and to evaluate (3 in (|2.31|) more 
explicitly. We begin with investigating the mass spectrum. We are interested in the modes whose mass eigenvalue 
m 2 is less than 2, since those modes contain a ghost (except for the case with m 2 exactly equal to 0). In the bulk 
Eq. (|2.18[) can be solved as 

Uj = B+y v + + B_y v - , (2.32) 

where B± are some constants and 



v± = 1/2 ± ^9/4 - m 2 . (2.33) 

From the conditions at y = y±, we obtain 



A 


B = 0, 


B := 






+ m 2 H. 


A: 


_((-» 






"U-f 


-m}H. 



I - m))HZ v+ -I - m)R_r c -J\ rn^Hl 



(2.34) 
(2.35) 



To obtain non-trivial solutions (B ^ 0), detyl = must be satisfied. Figure Q] shows the relation between iJ_ and m 2 
determined by det A = 0. Fig. Q] shows that no ghost mode appears if the relation 

H_> H s :=H + -r~\ (2.36) 

is satisfied. Otherwise, one ghost mode whose mass eigenvalue m 2 is smaller than 2 appears. 
Next we turn to the issue of finding a more explicit form for (3. For this purpose, we solve 



L (TT) + ^ J Q(y) + 4r c HlQ{y)5{y - y a ) = -JS(y - y+), (2.37) 

CT = ± 
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H 



FIG. 1: m as a function of H- 



in two different ways. First we solve this equation by expanding Q(y) in terms of Ui(y) as 

oc 



Substituting this expansion into the above equation (|2.37D and using the orthonormal conditions of the eigen- functions, 
we obtain 

Q(y) = ± u -^^J. 

mf — 2 



Another way is to solve Eq. (|2.37|) more explicitly. The bulk part is solved by 

Q(y) = C + Cia(y), 

and Cq and C\ are constants to be determined by the conditions at y = y± . These conditions are 



2H 4 



-1 + H+r c ) Co + (-1 + 2H+r c ) d = — - J, 



2iJ_ (1 + iJ_r c ) C* + (1 + 2H_r c ) C x = 0. 
Solving these for Co and Cx, we obtain 

-(l + 2F_r c ) 



(2.38) 



(2.39) 



C 

Ci = 



AH+(H + r c - 1)(1 + 2ff_r c ) - AH_{2H+r c - 1)(1 + H_r c ) 
H-(1 + H_r c ) 



J. 



J. 



(2.40) 



2H+(H+r c - 1)(1 + 2ff_r c ) - 2H_{2H+r c - 1)(1 + ff_r c 

Substituting (|2.40p into Eq. (|2.38|) , we obtain Q(y). Comparing the expressions for Q{y+) evaluated in these two 
different ways, we find 

1 2H_{1 + H_r c ) - H + (l + 2H_r c ) 



00 ii 



i=0 



mf - 2 H+ 4H + (H + r c - 1)(1 + 2tf_r c ) - AH_{2H + r c -!)(! + H_r c ) ' 



(2.41) 



Substituting this into Eq. (|2.31[) . we finally find 
3Hl(2H + r c - 1 



/3 = 



64K 2 (l + 2 J ff_r c ) 



(H+{H+r c - 1)(1 + 2ff_r c ) - H_{2H+r c - 1)(1 + #_r c )) 



3iJ2(2if +rc -i)2 H+(l + H. 



+ 'c) 



64k 2 



g_(l + g_r c ; 
(1 + 2ff+r c ) (1 + 2#_r c ) 



(2.42) 



where iJ + := iJ + — r^ 1 . One can see /3 changes its signature at H_ = H + . When H_ > H + (H- < H+), (3 is 
negative (positive). This means that only when the spin-2 ghost modes are absent, the spin-0 mode becomes a ghost 
mode. (See Eq. (|2.36[) .) Consequently, this model always has a ghost. 
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III. STABILIZATION 



In the preceding section, we have seen that when the spin-2 excitations have no ghost, the spin-0 brane-bending 
excitation becomes a ghost instead. Usually spin-2 excitations (i.e., transverse-traceless tensor modes) are completely 
decoupled from spin-0 excitations at linear order. However, we can construct a transverse-traceless tensor from a 
spin-0 mode whose mass squared is —AH 2 . For a scalar function $ that satisfies (D^ + 4)<f> = 0, we can compose a 
transverse-traceless tensor as 

«V - fv^V* - j7^ D ) *. (3.1) 

It is trivial to show that this tensor is traceless. The transverse property, = 0, can be shown by using 

Eq. (|2.17p . is also an eigen-function of cK 4 ), but its eigenvalue differs from that of $, which is —4. By using the 
identity (|2.26[) . the mass eigenvalue of the induced spin-2 perturbation for the operator (cK 4 ) — 2) is found to be 

m 2 = 2. (3.2) 

Usually, we identify H±m 2 as the mass of a rank-two tensor field seen by observers on the (±)-brane. Hence, the mass 
squared of this special mode is 2H± in the language of spin-2 perturbations. On the other hand, the mass squared of 
the corresponding spin-0 mode, which is the eigenvalue of H 2 0^ , is — 4H±, Here, we refer to the mass of this mode 
as the critical mass. This mass agrees with the threshold mass below which the spin-2 ghost appears. If the mass of 
a spin-2 mode is smaller than this critical mass, the helicity part of the mode becomes a ghost [f|. 

In the two-branes model discussed in the preceding section, the unique spin-0 excitation mode has exactly this 
critical mass. On the other hand, m 2 = 2 is the threshold spin-2 mass below which a ghost appears. If we gradually 
decrease the separation of the two branes, the ghost in spin-2 excitations eventually disappears. The mode that was 
previously carrying a ghost satisfies m 2 = 2 at the critical point. At this point this spin-2 mode can couple with 
the spin-0 excitations as an exceptional case. After passing the critical point, the ghost is transferred to the spin-0 
excitations. 

Then, what happens if we introduce a mechanism to stabilize the separation of the branes? Since the mass spectrum 
of spin-0 excitations changes (l9j |. the spin-0 mode with its mass squared being —AH 2 will no longer exist. In fact, 
the smallest mass squared must be large enough if the stabilization mechanism works. Then, again we consider the 
situation in which the brane separation is gradually decreased. When the spin-2 ghost disappears, we expect absence 
of the spin-0 mode at the critical mass, which can exceptionally couple to spin-2 excitations. Then there seems to be 
no reason why the spin-2 ghost is transferred to the spin-0 ghost at the critical point. 

However, by studying models stabilized by a bulk scalar field, we will find below that our naive expectation is 
wrong. What we will see is the following. With a bulk scalar field, there appears an infinite tower of Kaluza-Klein 
modes for spin-0 excitations. The smallest value of the mass squared of spin-0 excitations is not, in general, equal 
to —AH 2 as is expected. However, at the critical point where the smallest mass eigenvalue for spin-2 excitations m 2 
crosses 2, one of the values of the mass squared of spin-0 excitations also crosses the critical value —AH 2 . As a result, 
a spin-0 mode whose mass squared is smaller than —AH 2 arises if the smallest m 2 is greater than 2, and the spin-0 
mode becomes a ghost. We discuss how this happen in detail below. 

A. Background 

We introduce a bulk scalar field to stabilize the brane separation . The action is given by 

S = S 2 + S s , (3.3) 



S s 



er = ± 



By choosing the potential in the bulk Vb(V0 an< ^ the potentials on the branes V^(ip) appropriately, we can stabilize 
the brane separation. In the following discussion we do not need the explicit form of the potential functions. The 
unperturbed background configuration is similar to the previous case. The bulk metric is given by (|2.3[) and the 
branes are located at a fixed value of y(= y±) as before. But the functional form of the warp factor a(y) is different. 
The background bulk equations for the warp factor a(y) and the scalar field ip are 
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9V B n 



■4V , 
a oip 



(3.6) 
(3.7) 



Perturbations 



We use the "Newton gauge" , in which the spin-0 component of the shear of the hypersurface normal vector vanishes, 
following Ref. [19]. In this gauge, using the traceless part and {y/x}-component of the Einstein equations, we find that 
perturbations of the metric and the scalar field are related as 



Sip 



3 



2K 2 ip> 



[d y + 2H] (f>, 



(3.8) 



where 7 M „ := a 2 (y)^^ u , TL := a'/a and is a tensor which satisfies transverse-traceless conditions. Here, one 

remark is in order. For modes at the critical mass there is a possible coupling between spin-0 and spin-2 excitations. 
Recall that the traceless tensor generated by acting a derivative operator from a scalar-type function ([VpV„ — 
jj^n^'jffr ) is automatically transverse. Due to this mixing, vanishing of spin-0 component of the shear is slightly 
ambiguous. Here we define the "Newton gauge" by the above form of perturbations (|3.8[) . and therefore it actually 
can possess spin-0 component of the shear for the mode with this critical eigenvalue. In this sense, the gauge is not 
completely fixed. 



In the present gauge, the bulk equations for become 



iVT) + 1 (Q (4)_ 2) 



with h^ TT ^ as defined in (|2.15|) . While, the bulk equation for <f> becomes 

£(*) + -!_(□(*)+ 4)1 ^ = 0, 

with 



m := a 2 a, 



1 



d v a 2 



2k 2 



y a 2 ^ 2 y 



(3.9) 



(3.10) 



(3.11) 



We first write down the junction conditions to impose at the branes in the "Gaussian normal" coordinates. Per- 
turbation variables with an over-bar "~" are understood as those in the "Gaussian normal" coordinates, as before. 
The junction conditions for metric perturbations are 



±{dy-m)m 



T(±) 



1 



2k 2 

T -^-%vip'Sip + 2r c 



X„ 



1 



with X^y defined in (|2.10p . While, the condition for the bulk scalar field is 

±25$ = V" (±) {Tp)84>. 



(3.12) 



(3.13) 



The generators of the gauge transformation from the "Gaussian normal" coordinates to the Newton gauge take the 
form of 



tv - 

H±) ~ 



<t>{y')dy' = / MW + q ±) (x»), 



(3.14) 



y± 



r v (y'W / 4>Av"W' 



r v {y') 



dy'+t {±) (xP). 



(3.15) 



10 



Under this gauge transformation, the perturbation variables transform as 

h\t\y) = M») + 2V (M ^f ) - 2H% v (y)i 



v 

(±)' 



(3.16) 
(3.17) 



Substituting these relations into Eqs. (|3. 12[) and (13 . 1 3|) . we obtain the junction conditions in the Newtonian gauge. 
The conditions for the traceless part of metric perturbations become 



with H± := T~i(y±) and 



i(4) 



4 ) V 4 

On the other hand, the condition for the trace part of metric perturbations becomes 

^ (D (4) + 4)Z(±)=± ^ r(± ) ; 

where 

Z {±) = (lT2r c H±)i y {±) -r c 4>(y ± ). 
The junction condition for the scalar-type perturbation becomes 

2k 2 



(±)- 



± 



{Si> - V (1 T 2r.Hr 1 {Z { ±) + r c <t>)) = ^ (n (4) + 4) <f> 



at y = y± with 



yn{±) ±2ip"/ip' 



(3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



In the following discussion we assume > 0, which includes the simplest case of the tight binding limit; V"^ 



Combining the bulk equations and the junction conditions, we formally obtain the same equation as (|2.14|) for 
spin-2 excitations, and 



£<*> + _(□«>+ 4) cf>=J2 U^-(l-a2r c H ± )- 1 (Z (a) +r c ^--^(a^ + 4) ( f>U(y-y a ), (3.24) 



for spin-0 excitations. 



Solution with source 



As was done in £fl"fl we formally write down a solution of perturbation equations for a given energy momentum 
tensor by using the Green's function written in terms of eigen-functions of Z/ Tr ) and We also derive expressions 
for the metric perturbations induced on the branes. 

(TT) 

We begin with spin-2 excitations. We expand the metric perturbation h^v by using the eigen-functions of the 
eigenvalue problem (|2. 18|) . as before. The calculation closely parallels the previous case, although the functional form 
of the warp factor a is not the same. We obtain (|2.23[) . formally the same expression as before. 

Spin-0 excitations can be solved in the same way. We define the eigen-functions Vi{y) by 
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l+Y,2e a S(y~-y a ))v l (y) = 



(7=± 



—L 



(<t>) 



EAr c K 2 a 2 x 
cr (1 - o2r c H±) 6{y - y a 



and we define the inner product by 



dy 

^,2 



2e"5{y-y a ) \v i (y) Vj (y). 



Vi (y), (3.25) 



(3.26) 
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Under the assumption that e CT > 0, it is trivial to show that the norm of each eigenfunction is positive definite, and 
hence they can be normalized as (wj, Vj)w — Sij. Using thus normalized eigen- functions, the solution of Eq. (|3.24[) is 
given by 

it \ Vi(y+)vi(y) , , 

= ~ 2r c H+) D (4)_„2 Z (+)- ( 3 - 27 ) 

The induced metric on the brane for the solution with the source is given by 

*>$ = h% T) iv+) - %» (<KV+) + 2W+3V)) , (3.28) 
Here, transverse-traceless part h^F^ takes the same form as given in Eq. (|2.23[> . and it can be rewritten more explicitly 



as 



While the trace part is computed as 

2 / „,2( 



where for brevity we have introduced 



3Hl \3H* V Y M? + 4 y a) □ 



= tj2 (2rcH+ _ ^ _ (3 31) 



The above expression for hffl apparently contains a pole at the critical mass eigenvalue corresponding to (cK 4 ) +4) 1 . 
The term proportional to V M V„ in Eq. (|3.29|) is pure gauge. Hence, the contribution from the pole at the critical 

mass eigenvalue is proportional to 7 M „(n( 4 ) + 4)~ 1 T( + ). Collecting such terms in hffi coming both from the TT-part 
and from the scalar-part, the coefficient becomes 

In ijlll Dl we will prove an identity which shows that this coefficient vanishes as a whole. 



D. No physical degrees of freedom at critical mass 

As we mentioned earlier, our "Newton gauge" is not a complete gauge fixing. In this section, we will fist show that 
there is no physical mode at the critical mass in general after completely fixing the residual gauge. Next, we will show 
that the value of the expression (|3.32[) equals 0. 

The residual gauge degree of freedom of our "Newton gauge" exists at the critical mass, i.e., when cK 4 ) = —4 for 
scalar-type perturbations and cK 4 * = 4 for transverse-traceless tensor perturbations. To see this, let's consider the 
gauge transformation generated by the following infinitesimal coordinate transformation: 

e = ~a% (jj?) , e = -\r u x,»- (3.33) 

This gauge transformation keeps the conditions h yil — 0. The change of <fi induced by this gauge transformation can 
be read off the {y?/}-component of the metric perturbation: 



(3.34) 



12 



At the same time, the induced change of 4> can be also read off the trace part of the metric perturbation: 

H = -a'ad v (3-35) 
where we have used □WJC = — 4X . Equating these two expressions, we obtain an equation for X , which is 

1{TT) X = _ 2 x (3 36) 

a 2 

Gauge transformations with X satisfying this equation preserve the "Newton gauge" conditions. Hence, such gauge 
transformations represent the residual gauge degrees of freedom. 

Notice that the above equations are identical to the equations for the TT-part (|2 . 18[) with cK 4 ) = 4. The traceless 
part generated by this gauge transformation is 

Sh$P = - (V,V„ - |7^ (4) ) X. (3.37) 

Therefore, Eq. (|3.36|) means that the h^J" 1 induced by the residual gauge transformation satisfies the perturbation 
equation consistently. It is also straight forward to show that (j> = S(f> given by Eq. (|3.34[) , or equivalently by Eq. (I3.35|) . 
satisfies the equation for the scalar-type perturbation with D^ 4 ) = —4. To conclude, solutions of the bulk perturbation 
equations for both spin-2 and spin-0 parts at the critical mass are generated by a gauge transformation. 

Next, we turn to the the issue of boundary conditions. As we are interested in the gravitational degrees of freedom, 
we consider cases without matter fields: T^> = 0. From Eq. (j3~2U)) . wc find (D< 4 ) + 4)Z (±) = 0. Both spin-2 and spin-0 
parts are sourced by Z(±) and they have solutions for any given Z^±y However, the spin-2 part of the solution, which 
is proportional to (V M V„ — j^^n^^j Z(±y can be erased by the residual gauge transformation mentioned above. By 
the definition of residual gauge transformation, it does not violate the gauge conditions necessary to derive the set of 
perturbation equations ()2.14j) and l|3.24[) . Therefore the solution after this gauge transformation still satisfies all the 
perturbation equations. Once the TT-part is completely erased by the gauge transformation, the junction conditions 
therefore imply that Z(±j also vanishes after this gauge transformation. The equation of motion for the scalar part 
also continues to hold after the gauge transformation. But now, with vanishing -Z(±), the only consistent solution for 
(j) is cf> = 0, except for models with the potentials Vb an d being fine tuned. 

For any given homogeneous solution of -Z(±), we can construct a solution for the bulk perturbations which solves 
all the perturbation equations. However, what we have found is that seemingly non-trivial solutions constructed in 
such a manner are all pure gauge modes related to the residual gauge degrees of freedom of the "Newton gauge" . As 
we have anticipated, there is no physical perturbation mode at the critical mass. 

We found that by fixing the residual gauge, we can set Zt±\ = in the source-free case, and therefore we succeeded 
in eliminating the mode at the critical mass from the physical spectrum. Using this fact, we can prove that the 
expression presented in (|3.32|) vanishes as follows: By making a residual gauge transformation to the unperturbed 
background, we can set Zi + \ = Zq for such Zq that satisfies (cK 4 ) + 4)Zq = 0. After the gauge transformation, wc 

obtain a pure gauge mode for h/jfj"^ and (j> at the critical mass. From Eqs. (|3.35[) and (|3.37|) . this gauge mode is given 
explicitly by 

h^ = -{v^ v -\ lltu u^X i (3.38) 

<t>° = - a ' ad v (§)-§> ( 3 - 39 ) 

where a subscript c indicates the mode at the critical mass. Since it is a pure gauge mode, the above mode should 
satisfy the junction conditions (|3.22p and (|3.18[) with Z/ + -\ = Zq and T^t' = 0. Both conditions give the same 
boundary condition for X. Substituting (|3 .38[) into (|3 . 18[) at y = y + , we obtain 



a 2 d„ (4 



* a 2 



= 2Z Q -2r c H 2 + X\ y=y+ . (3.40) 

y=y+ 



Substituting this into Eq. (|3.39[) to eliminate d y X, we get 



MV+) - -2H+Z -Hl(l- 2r c H+) X. 



(3.41) 
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On the other hand, from Eqs. ((2~2"3")) and (j3~27f with the aid of ((3381) , we have 

1 = 4 Lr^ 4 and ^ = ^Ml^TTri z °- (3 ' 42) 

i 1 + i * 

Combining Eqs. (|3.41p and (|3.42p , we finally obtain 

Namely, the coefficient of the pole at the critical mass, (|3.32[) , vanishes. 

Although Eqs. (|3.29p and (|3.30|) apparently contain poles at the critical mass, we find that they completely cancel 
each other. This is consistent with the conclusion obtained earlier in this subsection that there is no physical mode 
at the critical mass. 



E. Effective action and ghost conditions 

Applying the same method used in Dl we now derive the four dimensional effective action to the quadratic order 
for the present case. From the derived effective action, we will identify the conditions for the presence of ghost. 

We denote the Kaluza-Klein towers of the transverse-traceless part and the trace part of the metric perturbation 
induced on the (+)-brane by h!^l and h!^\ respectively. The induced metric perturbation is given by 

i j 

Then, the kinetic term of the effective action is 

Skin = £>i [d i x^h^(n-2-m 2 i )h ( i }l 

i 

+ [ d 4 x^h^(a - M 2)fcCfl (3.44) 

and the action of the matter localized on the (+)-brane is 

S^ur = (0-th order of h) + £ ~ J dtxy/^foh® 

+ E I J d 4 xV=gT^h^ + (9(h 2 ). (3.45) 

j 

Comparing the equations of motion derived from the variation of Skin + S ma tter with Eqs. (|3.29p and (|3.30p . we 
find 



- 1 fl _ +4) 

" = ^uKy + y 128«V(» + )- 



All coefficients of the spin-2 component, a^, are always positive. However, when the mass eigenvalue is in the range 
< mf < 2, such a mode contains a ghost in its helicity zero component. In contrast, the coefficients of spin-0 
component, Pi, becomes negative when /if < —4. In this case this mode becomes a ghost. Hence, to realize a ghost- 
free model, all masses of spin-2 modes and spin-0 modes must satisfy mf > 2 and /if > —4. However, it turns out to 
be impossible that both conditions are satisfied simultaneously. 

In order to show this, we consider continuous deformation of the model potentials of the bulk scalar field. When the 
smallest eigenvalue of spin-2 excitations, m§, approaches 2 from below, UQ(y + )/(m§ — 2) in Eq. (|3.43p diverges. Since 
this equality is always satisfied, this divergence must be compensated by the other terms. The possibility is that one 
of the eigenvalues of spin-0 excitations /ij approaches —4 from above. Therefore, when the value of m\ exceeds 2, at 
least one of the values, /x 2 , must be smaller than —4. Hence, we cannot simultaneously satisfy m\ > 2 and /Xg > —4 
by a continuous deformation. All models specified by given bulk scalar field potentials are connected with each other 
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by continuous deformations. Thus it is proved that we cannot realize a model free from a ghost within the two-brane 
extension of the DGP model with a bulk scalar field. 

In the above discussion we have introduced the four dimensional effective action just from the notion of the equations 
of motion. Homogeneous equations of motion are sufficient to fix the form of the effective action except for the overall 
normalization of each mode. However, to determine if there is a ghost or not, the overall normalization is crucial . 
We have determined this overall normalization factor by investigating the coupling to the matter energy momentum 
tensor. 

There are alternative ways to fix this overall normalization for each mode without looking at the coupling to matter 
fields. We could have derived a four dimensional effective action directly from the original five dimensional action, 
although we feel that the approach adopted in this paper is slightly easier. One short cut to fix this normalization 
is just to look at the bulk part of the five dimensional reduced action written down in terms of the master variable. 
To obtain the five dimensional action, we can use double Wick rotation: y — > iy and t — ► —it, which transforms 
the four dimensional de Sitter space into an Euclidean four sphere. Then, the discussion in the bulk is completely 
parallel to the cosmological perturbation in the (five dimensional) closed FRW universe, and we can make use of the 
knowledge about the effective action in the context of the cosmological perturbation in the FRW universe [Hj]. Since 
our action should be obtained by the double Wick rotation from this action, one can easily read the correct form of 
the reduced action including the overall normalization. By doing so, we will find that when and only when [if < —4, 
the overall normalization factor becomes negative for spin-0 excitations. For spin-2 excitations, the overall factor 
is always positive definite. This gives an another route to reach the same conclusion that we have obtained in this 
section. 



IV. SUMMARY 

In the present paper, we have attempted to construct a ghost-free model by modifying the self-accelerating branch 
of the DGP brane world scenario, in which the spin-2 sector drives the accelerated expansion of the universe. We 
first tried adding an additional boundary brane in CTfl and then stabilizing the brane separation by introducing a 
bulk scalar field in mill In both cases, we found that the ghost excitation survives. As we have already mentioned in 
Introduction, this is caused by the degeneracy between spin-2 and spin-0 modes at the critical mass. This degeneracy 
enables the ghost property to transmute between these two normally decoupled modes. But this degeneracy occurs 
only at this special mass. 

In the case without stabilization, the fluctuation of brane separation gives the spin-0 mode. Hence, if we fix the 
brane separation, the mass spectrum of spin-0 excitations will be significantly modified. We introduced Goldberger- 
Wise mechanism [l8[ into our two-branes model. Then, in general the spin-0 spectrum does not have a mode at 
the critical mass. However, when the spin-2 mode crosses the critical mass, we found that the one of the mass of 
spin-0 mode also crosses the critical value, and it turns into a ghost. We have explained why we cannot avoid this 
exchange of a ghost degree of freedom in detail in the text. The point is as follows: If we compute the effective 
action for the matter field obtained after integrating out the gravitational degrees of freedom, the spin-2 contribution 
unexpectedly gives the term depending on the trace of the matter energy momentum tensor. This is because the 
spin-2 perturbation is sourced by the brane bending, which is sourced by the trace of the energy momentum tensor. 
Although this brane bending does not appear in the physical spectrum of the model once stabilization is introduced, 
its mass still corresponds to the critical mass. As a result, when one of the masses of the spin-2 excitations crosses 
the critical value, the pole of the brane bending resonates with the pole of this crossing mode. This leads to the 
divergence of the effective action of the matter field. Physically, however, any divergence in the effective action is not 
expected. This divergence is cancelled by the contribution from the spin-0 mode. But such cancellation is possible 
only if the mass spectrum of spin-0 modes also crosses the critical mass at the same time. 

It is instructive to emphasize the difference of the behaviour of spin-2 perturbations in the DGP model compared 
to the massive Fierz-Pauli model. In the massive Fierz-Pauli model, a pathology appears if the mass becomes 2H 2 
because the equation of motion implies that T = 0. In the DGP brane world, matter sources need not satisfy T — 0. 
Indeed, the solution for transverse-traceless perturbations shows no pathology at mf — 2. A pathological behaviour 
near m? = 2 in the spin-2 sector appears because we separated the spin-0 contribution from the spin-2 contribution. 
This is due to the degeneracy between spin-2 and spin-0 perturbations at m 2 = 2, which is explained at the beginning 
of mill This means that the DGP model has a natural mechanism to cure the pathology of the massive gravity 
theory at m 2 = 2 with T ^ by accommodating the physical spin-0 perturbation at m 2 — 2 when the mass of spin-2 
perturbations crosses m 2 — 2. However, we have shown that it is this mechanism that transfers the ghost between 
spin-2 and spin-0 modes. 

To conclude, we confirmed in this paper that it is really difficult to erase this ghost within the context of standard 
linear analysis. However, a different way to avoid appearance of ghost is proposed in Ref. [24 ] ■ We think still further 
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study is necessary to conclude whether this ghost is really harmful or not. 
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APPENDIX A: EFFECTIVE ACTION FOR SPIN-2 FIELD 



1. Exact Lagrangian 

To obtain the quadratic action for a massive spin-2 field, the method with the surest footing involves starting with 
the Einstcin-Hilbert action with cosmological constant and Fierz-Pauli mass terms. The Lagrangian up to quadratic 
order on de Sitter background is [6J 

L h = a^f-^nKxW* + iv^ A V^A + \ (Vft - 2V^H V„/i 

-\ (h^ + \h^ ^ (V^ - ^ 2 )) ■ (Al) 

Here the Hubble parameter is set to unity for simplicity. The overall normalization of the Lagrangian a can be 
different from the conventional one. (By redefinition of field, one can say that the coupling with the matter energy 
momentum tensor is unconventional.) The Lagrangian including coupling to the matter will be given by 

L = L h + L m , (A2) 

with 

L m = \^FlK v T» V . (A3) 

Equations of motion derived from the variation of this Lagrangian become 

/ m 2 \ 1 

a f A M „ - — ( hp, - J^h) J + ~T M „ = 0, (A4) 



with 



(V Q V^ a/3 - ah) - (V + ) ■ (A5) 



Since V^X^ = and VT^^ = 0, we have a constraint: 

vv = ( A6 ) 

Substituting back this relation to the trace of Eq. (|A4[) , we obtain 

h = -7T^4 — t^T. (A7) 
3a(m 2 - 2) v ' 
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Again, substituting the above constraint equations (|A6[) and (|A7[) into Eq. (|A4|) . we obtain 

a(n< 4 > - 2 - m 2 ) V=- T ^ - 3(^,2) (V M V V - (m 2 - 1) 7m „) T 

- (V - + (V„V, - ^ (4) ) T) - U Z[J:™\ ^T. (A8) 

In the last line, we divided the traceless part and the pure trace part. Comparing the solution of the above equations 
of motion with the result obtained from the original theory ()3.29j) . we can read the value of a. The contribution of 
the last term in (|A8[) is combined with the second term in Eq. (|3.29[) to obtain 

The first term in the square brackets is pure gauge, while the second term is cancelled by the pole at (cK 4 ) + 4) = 
in the trace part given in (13.30|) as was shown in £1111 Dl 



2. Simplified Lagrangian 

One may feel like to use the following simplified form of the effective Lagrangian with the constraints (|A6[) and 
2]) manifestly imposed by using Lagrange multipliers A M and A: 



This Lagrangian is not completely equivalent to as is shown below. Further disadvantage of this Lagrangian is the 
appearance of T^ v in the gravitational part of the Lagrangian. 
Equations of motion derived from the variation of L' h + L m are 

|(n( 4 > - 2 - m 2 ) V + \t»* - A (w „) + = 0, (All) 



1 

'■3a(mF- 2))~ 

1 



V"V = -V M o „ v „ 2 „^ T, (A12) 



Trace of Eq. (|ATT|) with the aid of Eq. (|A"l3)> gives 



6(m 2 - 2) 2 

While divergence of Eq. (|A11|) becomes 

□ (4) +4 _ TO 2 1 

T 

6(m 2 - 2) 2 

where we have used the equation 



□ (4) _ o - rn 2 1 

" ' T + -T - + 4A = 0. (A14) 



V p al Z, ^ T - - ( D< 4 )A M + V„V M A" ) + V M A = 0, (A15) 



VD (4 »/i F = (D^ + 5)V V - 2V^ 
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Eliminating V^A from the gradient of Eq. (|A14[) and Eq. (IA15[) . we obtain 
1 1 



2 (to 2 - 2) 

Divergence of the above equation becomes 

1 1 



V M (D (4) + 4)T - V^A" + 2dWa„ + 2V l/ V^A l/ = 0. (A17) 



D (4) ( D (4) + 4 ) T _ 3 ( D (4) + 4) y^A" = Qj ( A1 g) 



2 (to 2 - 2) 
where we have used 

V^n (4) A M = VV^V^" 

= (□(*) + 3)V Al A' 1 . (A19) 

If we neglect the possibility of adding an arbitrary solution of the homogeneous equation (Q( 4 ) +4)V Jll A^ = 0, we find 

□ (4) 

6(m 2 - 2) 

Substituting this equation into Eq. (|A14j) . we find 



= -ZTZ^ 7^T. (A20) 



_ 777 z — 1 

A = j—, -T. (A21) 

6(to 2 - 2) y ' 



Substitution of this equation and Eq. (|A20[) back into Eq. (|A15|) results in 



6(to 2 - 2) 



Again, if we neglect the homogeneous solution, we find 

A, = -^h^ V»T. (A23) 

If we substitute (|A"2"Tj) and (|A"2"3"1) into (jATTj) . we obatain the same equation as Eq. (|A8|) . However, strictly speaking, 
this system has propagating modes corresponding to the pole at (cK 4 ) +4) = 0, which we have neglected in obtaining 
Eqs. (jMT|) and (fA"2l|) . 

As is mentioned earlier, we need to include T^ v in the gravitational part of Lagrangian L' h . This is also unsatis- 
factory. If we remove from the Lagrangian L' h , we obtain the result with the pole at (C]( 4 ) + 4) = even if we 
neglect homogeneous solutions. 

This discrepancy between the simplified Lagrangian and the original one is not so strange. Here what we call 
constraints are not really constraints. In the Lagrangian formulation, they are just a part of equations of motion. In 
general, it is not always justified to impose a part of equations of motion before taking the variation of the action. 



3. remark on adding matter coupling to the reduced action 

In this paper, to derive the effective action with source term, we simplly added the universal coupling term 
(l/2)y/— 7 h^JT^ to the effective Lagrangian valid for the cases without source. Sometimes justification of this 
naive appoach is not so transperent for constrained system. In the context of cosmological perturbation, in order 
to handle constrained system, we often use the effective action written in terms of the physical degrees of freedom 
alone, which is obtained by imposing all the constraints on the action to eliminate unphysical degrees of freedom in 
the canonical foralism. In this case, even in the linear theory, it is not so manifest whether the resulting equations 
of motion derived from the action obtained by adding the coupling term (1/2)^/— 7 h [lv T llv to the already reduced 
Lagrangian are identical to those obtained by taking the coupling term into account from the beginning. Here we 
give a proof of this identity. 

First, we briefly describe the method for deriving the reduced action in the source-free case. We consider the 
quadratic Lagrangian L(qi,qj) where qi represent the dynamical variables. The corresponding Hamiltonian and the 
primary constraints are, respectively, denoted by H{qi,pj) :— J^jPjQj ~ -k(<7i;<?j) and Cj = J2j a j ( lj + J2k PkPk = 0, 
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where pj :— dL/dqj. a,j and 0k are, in general, given functions determined by the background, but we here assume 
that they are just constants for simplicity. Following the Dirac's prescription [? ], we make a "new" Hamiltonian by 
adding primary constraints with Lagrange multipliers A^ as 

H> = H(Qi, Cj) + Wt, ( A24 ) 

i 

where Qi{qi,Pi) are the combinations of variables which are independent of CV From the consistency conditions for 
the time evolution of constraints, the secondary constraints follow. We denote them by Xi- Here we assume that 
the set of constrains C\ including both C and Xj is second class. Namely the matrix composed of the commutators 
between constrains is non-degenerate 

det([Cj,C$])^0. (A25) 

If the set of constraints is not second class, we add some gauge fixing conditions to make it second class. Notice that, 
from the consistency conditions, we also have equations for A^. Hence, now Xi are given as functions of variables Q[, 
where Q\ are the combinations of variables which are independent of Cj. (One can add Cj to A^ without destoying 
consistency conditions.) Thanks to (|A25[) . we can define Q\ in such a way that 

[QJ,C$]=0. (A26) 

Here a key point is as follows. To obtain a consistent set of constraints, in general, we need to add new term to 
the Hamiltonian corresponding to the secondary constraints Xi as we did for the primary constraints in Eq. (|A24[) . 
However, in practice, these additional terms are not necessary in most cases that we are interested in, although general 
proof is lacking as far as we know [25| . Here we assume this property. 

We rewrite the Hamiltonian as 

H' = H(Q' t , Cj, Xk) + E MQi)C*. (A27) 

i 

The consistency conditions for the constraints mean that [C,-, H'j does not contain Q' rj . This fact with the aid of (|A25[) 
and (|A26p means that H' does not contain terms in the form Q'fi'y Since AjCj has only terms in the form of Q^Cj, 
we find that H does not have terms in the form of Q^Xj- The reduced Lagrangian is obtained by setting C[ = in 

Y,Pi{Q'v C kMQ'vC' k ) - H'{Q' v C' k ). (A28) 

i 

Here the first kinetic term also does not have any cross-term between Qi and Xj because Qi are chosen so that they 
satisfy Eq. (|A"2"6]) . 

We are ready to proceed to the setup in which the source term, J^. qiSi, is present. To distinguish the quantities 
defined under the presence of the source term, we associate an overbar " ~ " with them. Lagrangian becomes L = 
L(qi, qj) + J2k QkSk- Adding source term ^ qiSi does not affect the definition of the conjugate momentum and hence 
that of the primary constraints. Hence they are identical to the source-free case, 

Pi := dl/d qi = pi, Cj = Cj. (A29) 

The modification of the Hamiltonian is also trivial. The "new" Hamiltonian with constraint terms becomes 

H' = H (<& C,-,Xk) - Yl 1 (Se + E ^ ( A3 °) 
i i 

As in the source-free case, secondary constraints "xi" follow from the consistency conditions. Since the limit S% — > 
should recover the source-free case, in linear theory the secondary constraints must take the form 



Xi 



Xi + ^ctijSj, (A31) 



where otij are some constants. For the same reason the Lagrange multiplier "Aj" also become A^ = A,; + J2jPijSj- 
The reduced Lagrangian is obtained by setting C' = in 



^PniQ'i, Cj,Xk(xe, S m ))q n (Q'i, Cj,Xk(xt, S m )) 

n 
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-H (Q'i, Cj, X k(Xi, S m )) + InSn + \(QIS 3 )C P . (A32) 

n p 

As we have shown earlier, the first kinetic term and H have no cross-term between Q\ and Xj- Therfore, these terms 
contain no cross-term between Q[ and Si, either. In the above expression for the Lagrangian, only the term q n S n 
contains the coupling between Q[ and This means that when we consider the system with the source term, it is 
justified to add the source term to the reduced Lagrangian obtained for the source-free system. 

A little subtlety arises when we reconstruct all the perturbation variables from the physical degrees of freedom 
Q'y To do the reconstruction, we need to solve the constraint equations. Here we need to recall that we are using 
a different set of constraints in two cases with and without source. As was mentioned above, Xi is n °t identical to 
Xi- In this sense therefore the difference between these two cases is not completely captured by the difference in the 
reduced Lagrangian. However, what we do in this process is just to solve equations which does not have temporal 
differentiation. Therefore even if we use x% = in error instead of Xi = in the case with source, the error is such 
that can be absorbed by a shift of qi to qt + Y) j 7,-j Sj . 
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